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Abstract 

The consequences of taking the generalized Chaplygin gas as the dark energy constituent of the 
Universe on the gravitational waves are studied and the spectrum obtained from this model, for the 
fiat case, is analyzed. Besides its importance for the study of the primordial Universe, the gravi- 
tational waves represent an additional perspective (besides the CMB temperature and polarization 
anisotropics) to evaluate the consistence of the different dark energy models and establish better 
constraints to their parameters. The analysis presented here takes this fact into consideration to 
open one more perspective of verification of the generalized Chapligin gas model applicability. Nine 
particular cases are compared: one where no dark energy is present; two that simulate the A-CDM 
model; two where the gas acts like the traditional Chaplygin gas; and four where the dark energy 
is the generalized Chaplygin gas. The different spectra permit to distinguish the A-CDM and the 
Chaplygin gas scenarios. 
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1 Introduction 

The most recent results from type la supernovae observations [I] |3 OJ and the cosmic microwave back- 
ground anisotropics detection ^ have led cosmologists and astrophysicists to the conclusion that most of 
the matter of the Universe interacts in a repulsive manner (as an example through a negative pressure). 
There are many candidates to describe this exotic fluid. The most natural candidate is a cosmological 
constant, but it presents a discrepancy of 120 order of magnitude between the theoretical predictions and 
the observational data 5 j. A self-interacting scalar field, known as quintessence, is another proposal to 
explain dark energy jSj. However, it asks for fine tuning of microphysical parameters in order to have 
a suitable potential term 0. Among many other possibilities, the Chaplygin gas model has receveid 
recently special attention IHUHinSl- The Chaplygin gas model is based on a perfect fluid whose pressure, 
besides to be negative, varies inversely with density. One of the interesting aspects of this fluid is a con- 
nection with branes in the context of string theories . Some phenomenological generalizations 
of this fluid have been proposed, leading to the so-called generalized Chaplygin gas model j^j. 

The aim of this paper is to investigate the particular signatures of the generalized Chaplygin gas, 
mainly in comparison to the standard A-CDM, in what concerns gravitational waves. Special attention 
will be given to the spectral distribution of energy density. 
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The gravitational waves are very important to Cosmology. They must have special signatures in 
the polarization [H] of the CMB anisotropics [T^j. Moreover, since gravitational waves has decoupled 
from matter already in the deep early Universe, they can be a window to the primordial phase. Even if 
gravitational waves (and consequently the polarization of the CMB photons x ) have not been detected 
directly until now, great efforts are been done in this sense and there is a hope that the next generation 
of experiment in space (LISA, Planck, etc.), may allow this detection in perhaps ten years. The GW 
spectrum frequency range of observational interest extends from 10 -18 Hz to 10 10 Hz and the energy 
density spectrum is constrained by the CMB, in units of the critical density &gw, as QH 
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This constraint will be used later in the evaluation of the spectra. Many works have been made in order 
to identify specific signatures of cosmological models in the spectra of gravitational waves, for example, 
in the case of quintessence model [TBI and string cosmology [H3 EH 12] • 
The generalized Chaplygin gas is characterized by the equation of state 

P-~ • <»> 

where A and a are constants such that A > and < a < 1. If the energy- momentum tensor conservation 
is taken into account, the relation between the generalized Chaplygin gas density and the scale factor a 
becomes ( 

where B is an arbitrary integration constant. The main properties of the relation (|HJ) are: (i) it interpolates 
the cosmological constant phase (when the scale factor is large) and a pressureless fluid (for small values 
of a) phase; (ii) the sound velocity associated with it goes from zero, for A = 0, to the velocity of light, 
being always positive. 

We will analyze the gravitational waves spectrum obtained from the generalized Chaplygin gas model, 
which is described in section O In section [3| we discuss general properties of the gravitational waves 
differential equation in the context of two important dark matter models: cosmological constant and 
Chaplygin gas. The spectrum obtained from numerical calculations on these cases are presented in 
section 0] and analyzed in section where we also present our conclusions. 



2 Outline of the model 

We consider a flat, homogeneous and isotropic Universe described by the Friedman- Robertson- Walker 
metric, which may be written as 

ds 2 = c 2 dt 2 - a 2 {t) gff dx i dx j , gff = ^ , (4) 
and leads the Einstein's equations to assume the form 
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where a is the scale factor of the Universe, while p m and p c are the pressureless fluid and the Chaplygin 
gas densities. The pressures p m and p c of the fluids are related to their densities by the equations of state 
p m = and p c = -A/p%. 

1 There are two polarizations modes, called E and B. The mode E has already been identified, but only the detection of 
the mode B will allow to identify gravitational waves 1161 1171 . 
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The option of working in a flat Universe is justified by the recent data from the CMB measurements. 
The curvature of the Universe is characterized by the parameter Qk which is defined, in terms of the total 
observed density pr and the critical density p cr , as Vtk — 1 — ^t, = Pt / per- The CMB data gives 

si k = o±o.o6 E2J. 

If the fluids interacts only through geometry, the energy-momentum tensor for each component is 
conserved and we get 



Pm = ^ , Pc = ( A 



B 



a 3(a+l) 



(7) 



We take the scale factor today as the unity, ao = 1, and thus p mo and p Co = (A + B) "+ 1 are the present 
densities of the fluids (The subscripts o, according to the current notation, indicate we are considering 
the present values of these quantities). From this last equation we can express the integration constant 
B in terms of A, B = p^ 1 - A. The Chaply gin gas density is then rewritten as 
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and the parameter A is connected to the sound velocity in the gas, v s , by the expression 
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dp c 
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Taking the first of the equations J7J) ( which refers to the pressureless fluid ) and the equation (|HJ), we 
use the relation Q io = pi /p cr (where i = m, c) in J5J) to obtain 
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where the Hubble constant Ho is defined by the expression Ho = do/do- Since we are restricted to a 
flat Universe, the fractions of pressureless matter and Chaplygin gas today, ^ mo and ^ Co , obey to the 



relation S^L 



i. 



With these two last equations we are able to write the GW amplitude differential equation as a 
function of the observable variables i^o, Q mo and ^ Co , and of the Chaplygin gas parameters A, a. It's 
also important to remark that if A = 0, the gas behaves like the pressureless fluid (and the situation is 
the same as if we had set ^ mo = 1) while, on the other hand, it behaves like the cosmological constant 
fluid when A = 1 (and therefore we can simulate the A-CDM scenario). 

Among the many possible cases produced by the combinations of parameters we chose a few important 
ones and classify them according to the fractions of dark energy and matter and to the kind of dark energy, 
as shown in tabled below. 

In the first case no dark energy is present, while the other eight cases reflect the situation where the 
dark energy constitutes all the non- visible matter (that means, no dark matter exists and ^ mo = 0.04) 
or the one where the dark matter and the visible matter altogether represent 30% of the total energy 
density (^ mo = 0.3). 

After establishing the characteristics of the studied models, we start in section|31the description of the 
behavior of gravitational waves due to different cosmic fluids contents, namely the cosmological constant 
and the generalized Chaplygin gas. 
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Table 1: Models of interest for gravitational waves power spectrum analysis. 

3 GW equation in Chaplygin gas and A-CDM models 



Cosmological gravitational waves are obtained by means of a small correction hij on equation Q. which 
represents the metric. Hence, the tensor related to the unperturbed metric, is replaced by gij = 
gff + and the resulting expression is [2H EU • 
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where /c is the wave number times the velocity of light (k = 27rc/A), the dots indicate time derivatives 
and we have written, hij(t,x) = h(t)Qij, where Qij are the eigenmodes of the Laplacian operator such 
that Qa = Qki,k = 0. 

Performing a variable transformation, from time to the scale factor a, and representing the derivatives 
with respect to a by primes, equation ([T2|) assumes the form 
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By using the background equations (fTUj) and (fTT|). or similar ones concerning to other fluids, into JT3J), 
one can easily express ft in terms of the parameters of the Chaplygin gas or any other model. 



3.1 Chaplygin Gas 

Let us perform the last operation mentioned above and find ft as a function of the redshift z, with 
the following steps: (i) use the fairly known relations l + z=^-,ao = l; (ii) perform a second variable 
changing (from a to z)\ and (iii) take back the dots to indicate, from now on, the new integration variable. 
These operations result in 
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where k has been redefined to absorb the Hubble constant, i.e., k = 2ttcHo/X. Therefore, h is a dimen- 
sionless quantity. 

The solutions h(z) of the differential equation oscillates with greater amplitudes as z — > 0: these 
amplitudes are increasing from the decoupling era up to the present. 

Equations (|14I16|) above may be used to perform all the cases of interest mentioned before. Setting 
Q mo = 1 and Q CQ = 0, for example, we have the pressureless fluid equation and for the cosmological 
constant we set A = 1. 



3.2 Cosmological constant 

In order to establish the GW equation for the cosmological constant fluid, we use the fact that the 
Chaplygin gas reproduces its behavior if A = 1, for any value of a, as stated before. Hence, the equations 
JTHJ) and (JTfljl become 



fl — ^71 



(17) 
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The resulting equation is easier to integrate than the one obtained before and presents the same 
qualitative behavior with respect to the redshift. Another important feature to be mentioned is the 
existence of quite simple analytical solutions, which permit us write h(z) in terms of a combination of 
sines and cosines: 



h(z) = 
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k 3 / 2 (l + z) 3 



(k(l + z)Cx + C 2 ) cos [k(l + z)\ + (k(l + z)C 2 - Ci) sin [fc(l + z)\ 



(19) 



The arbitrary constants C\ and C2 remain undetermined since the set of initial conditions is unspec- 
ified. Nevertheless, it is important to remark that, for a fixed /c, the argument of the cosine and sine 
functions in JTJJ) is = k(l + z) and so the oscillation frequency decreases when z tends to 0. 

Figure {Q illustrates the case where h(zi) = h(zi) = 10 -6 , z% = 4000 are the initial conditions for two 
different values of k. 




Fi gure 1: Gravitational waves curves for A-CDM model. The dash line refers to k = 1/100 and the solid line to k = 1/50. 



The behavior described is evident and one can also observe that the amplitude grows with time (i.e., 
as z — » 0), as we have remarked for the Chapligyn gas. This result implies that initial small primordial 
fluctuations on the gravitational field are amplified. 
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4 GW spectra 



The power spectrum of gravitational waves, defined as [T8] 
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(where ho(v) = h(0) and v = H§k /2ir), is generally obtained directly from the solutions of JT3J). In the 
particular case of the A-CDM model, an analytical result is also possible, as shown in and this fact 
may be used to verify the accuracy of the calculation and the applicability of the algorithm. 

For each of the cases of interest found in table QJ, we have assigned some common parameters, namely 
the initial conditions h(zi) = ^ -3 10 -5 , Zi = 4000; the range of frequencies 10 -18 Hz < v < 10 _15 i7z, 
and the normalization constant imposed by the constraint equation QJ. The factor v~ 3 in the initial 
conditions fix the primordial spectrum. The resulting spectra are presented in figures (|2I12|) . 
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Figure 2: Graphic (dn G w I dlnv) X {y X 10 _18 ifz) for 



the model where fl c 



and Or, 



Figure (|2j) refers to the pure pressureless fluid case and one can see that for v — > the amplitude 
reaches its maximum. It decreases very quickly and, for greater values of ^, it increases slowly. 

On the other hand, the spectrum corresponding to a cosmological constant dominated Universe JSJ 
has a different behavior. The amplitude grows faster as the frequency increases and reaches a value six 
times greater than the previous one. Figures © and (0 shows the behaviour of the spectrum for very 
low frequencies for the A-CDM case. For the other models, the spectrum is very similar. Still referring 
to the A-CDM model, figure (@J) shows a behavior similar to the one from the pressureless fluid, but the 
amplitude growing is faster. 
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Figure 3: Graphic (dn GW /dlnis) X {y X lO" 18 ^) for Figure 4: Graphic (dn GW /dlnis) X {y x lCr 18 /fz) for 



the model where A = 1, fl c 



-- 0.96 and ttn 
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■■ 0.7 and Qn 



■■ 0.3 



6 



4X10" 



£x!0 



1x10 




5x10 
5X10' 
4X10' 
3x10" 
£xlo' 
1X10' 




10 £0 30 40 50 



10 £0 30 40 50 50 



Figure 5: Graphic (dVL GW /dlnv) X (y X 1CT 18 Hz) for 
the model where A = 1, Q Co = 0.96 and Q mo = 0.04 for 
small values of k. 



Figure 6: Graphic (dQ GW /dlnv) X (i/ X KT 18 ifz) for 



small values of fc. 
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Figure 7: Graphic (dQ GW /dlnv) X (ux Kr 18 if z) for the 
model where A = 0.5, a = 1, £7 Co = 0.96 and £7 mo = 0.04 



Figure 8: Graphic (dn GW /dlnv) X {y X 10- 18 ifz) for 
the model where A = 0.5, a = 1, r2 CQ = 0.7 and Qm = 0.3 
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Figure 9: Graphic (dQ GW /dlnis) X {y X lO" 18 ^) for Figure 10: Graphic (dQ GW /dlnis) X {y x lO" 18 ^) for 



0.04 



0.3 



The other graphics correspond to the Chaplygin gas and present a similar shape. However, in general 
the Chaplygin gas models lead to a smaller amplitude in the gravitational wave spectra. This amplitude 
is still smaller for the case of pure dark matter model. Its also important to remark that figures (|7I8|) 
refers to the Chaplygin gas in its traditional form (with the equation of state p = —A/p, since a = 1) 
while the others correspond to the generalized Chapligyn gas. 
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Figure 11: Graphic (dn GW /dlniy) X (y X lO" 18 Hz) for 
the model where A = 0.5, a = 0, f2 Co = 0.96 and r2 mo 
0.0 1 
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Figure 12: Graphic (dn GW /dlnis) X {y X lO" 18 ^) for 
the model where A = 0.5, en = 0, fl CQ = 0.7 and flm = 0.3 



5 Conclusions 

In this work we have studied the fate of gravitational waves in the context of cosmological models where 
dark energy is described by the generalized Chaplygin gas. We have exploited the cases where dark 
energy and dark matter are unified by the Chaplygin gas, as well as the cases where a pressureless dark 
component is present besides the Chaplygin gas. In all situations, however, we have taken into account 
the baryonic component as infered from the primordial nucleosynthesis. We perform also an analysis, 
at same time analytical and numerical, of the A-CDM model and of a pure dark matter model (where 
pressureless matter is the only component of the Universe). In this sense, the goal of the present work 
was to try to identify special signature of each model in the energy spectrum as function of the frequency. 

The spectra for all models exhibit the same shape: initially the energy decreases with frequency, 
becomes almost constant, and then it increases with frequency. However, from the results it comes out 
that the A-CDM model presents a spectra with an amplitude greater then the pure dark matter model by 
a factor of 6. The Chaplygin gas model interpolates these two cases. In a typical situation, for example 
when a = 1, A = 0.5 and there are only Chaplygin bas and baryons, the amplitude of the spectra is about 
4 times smaller than in the A-CDM case. As the A-CDM case is approached (A — > 1), the amplitude 
increases, but in a very low rate, exceptly at the neighborhood of A = 1. 

These results show that gravitational waves discriminate very poorly the different models. Somehow, 
this is natural since the gravitational waves equation is sensitive essentially to the background behavior, 
in contrast with density perturbations which depends strongly on the kind of matter content. This is 
more accentuated in the comparison of different generalized Chaplygin gas models, that is, models with 
different a. In fact, the Chaplygin gas models are more sensitive to the parameter A then a. 

However, it must be stressed that even if the different models based on the generalized Chaplygin 
gas may be quite degenerate in what concerns the behaviour of gravitational waves, it is possible to 
distinguish them from the A-CDM model by analyzing the typical amplitude of the spectra. But, a 
complete study of the value of this amplitude asks for a complementation of the analysis here done 
by introducing the primordial spectrum of gravitational waves, for example, that one coming from the 
primordial inflationary phase, which allows to fix preciselly the amplitude and the initial power spectrum 
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